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Abstract 

The integrative and conjugative element ICEc/c Is a mobile genetic element In Pseudomonas knackmussii B13, and an 
experimental model for a widely distributed group of elements in Proteobacteria. ICEc/c Is transferred from specialized 
transfer competent cells, which arise at a frequency of 3-5% In a population at stationary phase. Very little Is known about 
the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene 
operon encoded by ICEc/c, which exerts global control on transfer Initiation. The operon consists of three consecutive 
regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that 
MfsR autoregulates expression of the operon, whereas TclR Is a global activator of ICEc/c gene expression, but no clear role 
was yet found for MarR. Deletion of mfsR Increases expression of tciR and marR, causing the proportion of transfer 
competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two 
orders of magnitude loss In population viability, Individual cell growth arrest and loss of ICEc/c. This indicates that 
autoregulatlon is an Important feature maintaining ICE transfer but avoiding fitness loss. Blolnformatic analysis showed that 
the mfsR-marR-tciR operon is unique for ICEc/c and a few highly related ICE, whereas tciR orthologues occur more widely in a 
large variety of suspected ICE among Proteobacteria. 
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Introduction 

Comparisons between ever-increasing numbers of sequenced 
genomes reveal the large extent to which prokaryotic genomes 
have undergone horizontal gene transfer (HGT) [1-5]. HGT has 
traditionally been viewed as the consequence of natural transfor- 
mation, or of the action of mobile elements such as conjugative 
plasmids and phages [6,7]. During the last decade, however, other 
types of mobile genetic elements such as integrative and conjugative 
elements (ICEs) have been recognized, which are widespread and 
thus may significantly contribute to HGT [8-13]. In contrast to 
phages and plasmids, however, we stiU know little about the life 
styles of the diverse ICE types, their modes of self-transfer and 
regulatory pathways controlling self-transfer. Like temperate 
phages, ICEs mostly exist in an integrated form at one or more 
specific sites in the host's chromosome (often in genes for tRNA), 
and are vertically transmitted to daughter cells by chromosome 
replication and segregation [8,10,14]. In order to transfer horizon- 
tally, ICEs excise themselves by site-specific recombination {attL and 
attR, Figure 1). This produces a circular double-stranded DNA 
molecule, which can transfer by conjugation to a new recipient cell. 



where it can reintegrate [14]. Autonomous plasmid-like replication 
of the excised form may occur [15-17], but is not required for the 
transfer itself. 

The regulatory mechanisms that control the switch from 
integrated to excised state vary widely among different ICE 
types insofar as this has been studied. In several ICEs, this 
switch is the consequence of a cascade of a variety of regulatory 
factors, such as Phrl/Rapl and ImmR/ImmA in ICE&7 [18], 
SetR/SetCD in ICESXT [19,20], KorSA/Pra in pSAM2 [21- 
23] or QseM/TraR in ICEMlSym^'^ [24,25]. Most wild-type 
ICEs transfer at low frequencies (i.e., less than 1 per 10'^ donors), 
suggesting that the regulatory cascades keep extremely tight 
control and allow only a small subset of cells in a population to 
follow the path of ICE excision and transfer, but the need for 
such tight control is a priori unclear. This bistability is most 
pronounced and well-studied for a model ICE named ICEc/t in 
Pseudomonas [26,27], which is evolutionary very distinct from the 
afore-mentioned ICEs [8,28]. ICEc/c is originally found in two 
copies in Pseudomonas knackmussii B13 and is member of a family 
of ICEc/c-like elements widely distributed among proteobacterial 
species [29]. ICEc/c is integrated at the 3'-eiid of tRNA'''-' genes 
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Author Summary 

Integrative and conjugative elements (ICEs) are a rela- 
tively newly recognized class of mobile elements in 
bacteria, which integrate at one or more positions in a 
host chromosome, can be excised, circularized, and 
transfer by conjugation to a new recipient cell. Genome 
sequencing indicated that ICEs often carry genes with 
potential adaptive functions for the host. Various ICE- 
types have been described and ICEc/c is a useful model 
for a wide class of elements found in Beta- and 
Cammaproteobacteria. Because ICEs normally remain 
"silent" in the host chromosome and often lack selectable 
markers, their lifestyle is difficult to study. One of the 
characteristics of ICEc/c is that transfer is initiated in only 
3-5% of donor cells in a population during stationary 
phase. Here, we describe an operon of three regulatory 
genes, two of which control the transfer initiation of 
ICEc/c. Our findings suggest that the low transfer rate 
results from the repression of an activator and that this is 
essential to minimize the deleterious effect of hyper- 
activation of transfer initiation. While the individual 
regulatory genes are quite common on ICEs, they rarely 
occur in this configuration. 



but can excise itself by the action of the IntBlS integrase 
encoded on the element (Figure lA). Expression of intBlS in the 
integrated form is under control of the promoter Pjji,, which by 
single cell reporter gene analysis was shown to become active 
only in 3-5% of a bacterial population during stationary phase 
[26]. Direct single cell visualization further confirmed that only 
cells which express reporter gene fused to Pi,,; above a threshold 
are capable of transferring ICEc/c to new recipients, a bistable 
state which we recently named "transfer competence" (tc) [30]. 
Irrespective of the success of ICEc/c transfer, tc cells can only 
divide a few times once they re-enter exponential phase before 
they arrest growth. We recently showed that this is due to the 
expression of the ICEc/c genes shi and parA [30] . Expression of 
intBlS is dependent on a variety of factors, most notably a gene 
named inrR (Figure lA), which itself is also bistably expressed 
[26]. Both inrR and intB13 expression are dependent on the 
abundance of the stationary phase sigma factor RpoS, with cells 
having highest RpoS levels being more likely to activate P;,,, and 
PinR [29]. RpoS and InrR are important for activating ICEc/c 
excision and transfer, but are not sufficient. Therefore, we 
hypothesized that additional factors are necessary for the tc state 
to develop [29]. 

In this study, we report a locus of three consecutive regulatory 
genes on ICEc/c, which is essential for controlling its transfer. 
The locus was uncovered by random transposon mutagenesis, 
and further studied by creation of deletion mutants and 
complementation. The effect of mutations was studied at the 
level of ICEc/c expression through microarray hybridizations, 
RT-PCR and reporter gene-based single cell fluorescence 
microscopy, and further in ICEc/c transfer assays. Fitness of 
mutants compared to wild-type was examined in growth assays 
and individual cell fates were followed by microscopy. Bioinfor- 
matics was used to analyze the configuration of the ICEc/c 
regulatory locus within this ICE family, and to possibly 
reconstruct the steps that may have led to selection of the 
specific regulatory control mechanism of ICEc/c. The results of 
our study help to explain why a careful balance has to be 
maintained between ICE transfer frequency and fitness loss. 
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Figure 1. Schematic overview of ICEc/c and the location of the 
genes relevant to this study. (A) ICEc/c integrated form with the two 
flanking 18-bp repeats (as black rectangles, attR and attL). Previously 
determined gene regions are placed at their approximate location. 
Important functional regions are depicted as grey boxes accompanied 
by legends. + or -, indicate the orientation of the coding strand for the 
respective gene(s) (the + strand being the one of intB13). dc genes, 
chlorocatechol degradation; amn genes, 2-aminophenol degradation, 
kb, kilo-basepairs. (B) Detail of the mfsR operon. Arrows, predicted open 
reading frames (the right-to-left orientation indicates the minus strand). 
Triangles, positions of the Tn5-mediated kanamycin gene insertions 
(nucleotide positions indicated below, according to the AJ61 7740.2 
numbering). Regions deleted in this study are displayed as white bars 
with the names of the mutations noted below. The hooked arrow 
indicates repression by the mfsR gene product on the P^fsp promoter. 
(C) Detail of the fragments inserted by mini-Tn7 delivery for 
complementation of mfsR and tciR deletions. (D) BlastP-predicted 
domains for each of the three regulatory genes in the mfsR operon. C 
and N, carboxy and amino terminus, respectively. 
doi:1 0.1 371/journal.pgen.l 004441 .gOOl 



Results 

Discovery of an ICEc/c transfer control locus by 
transposon mutagenesis 

In order to discover ICEc/c-located factors involved in its self- 
transfer, a library of P. knackmussii B 1 3 mutants was generated by 
using random Tn5 mutagenesis [31]. Next, we recovered ICEc/c 
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elements with Rm-insertions by conjugating the pool of B13 
mutants en masse to Pseudomonas putida UWC 1 and selecting for Km- 
resistant P. putida (Figure SI). We hypothesized that mutant ICEcfc 
with insertions in genes implicated in self-transfer could still be 
transferred to UWCl, when the second copy of ICEc/c in the same 
B 1 3 donor cell is intact and complements transfer of the mutant 
copy. A total of 1920 Km-resistant P. putida transconjugants was 
recovered and subsequendy conjugated each individually with a 
second P. putida recipient, resistant to nalidixic acid (Figure SI). 
For those conjugations in which no Km- and nalidixic acid- 
resistant transconjugant growth was detected, the corresponding P. 
putida donor was recovered and the location of the Km'^^-gene 
insertion on ICEc/c was mapped. A total of 18 clones was 
recovered, which had insertions in an ICEc/c open reading frame 
numbered orfl8502, that we renamed mfsR (Figure IB). Surpris- 
ingly, apart from one donor with an insertion in intB13, no other 
mutants with impaired ICEc/c transfer were found in this 
screening. The Km^-gene had been inserted in four different 
positions in mfsR, at ICEc/c nucleotide positions 19033, 18758, 
18730 and 18618 (Figure IB, accession number: AJ61 7440.2). 
This suggests that transposon insertions in strain B13 were 
suSiciently frequent to cover all genes, but that the selection 
procedure was biased for the recovery of the mfsR insertion, which 
may have been due to the function of mfsR as regulator in ICEc/c 
transfer (see below). Alternatively, it is possible that insertions 
in ICEc/c genes needed for transfer might not be efficientiy 
complemented by the second ICEc/c copy and would thus be 
underrepresented in the P. putida library. Frequencies of ICEc/c 
transfer of the strains P. putida UWCl -2961 (Km'^-gene insertion 
at 19033) and UWCl-2962 (insertion at 18618) in a filter-based 
conjugation assay were lO'^-fold and 10^-fold lower than of a P. 
putida with one integrated wild-type ICEc/c copy, respectively 
(Figure 2A). 

mfsR is part of an operon formed by three consecutive 
transcriptional regulators 

Closer inspection indicated mfsR to be the first open reading 
frame in a series of three consecutive transcriptional regulators, 
previously designated as orfl8502, orfl7984 and orfl7162 
(Figure IB). mfsR encodes a TetR-like regulator harboring helix- 
turn-helix motifs TetR_N and TetR_C_7 (pfam0040 and 
pfaml4246, respectively, see Figure ID). The orfl7984 gene 
overlaps with the end of the mfsR open reading frame by 4 bp and 
encodes a putative regulator of the MarR family (smart00347 
HTH_MARR motil). The last gene of this cluster starts 24 bp 
downstream of the stop codon oi orfl7984 and is predicted to code 
for a LysR-type transcriptional regulator, harboring an N-terminal 
HTH_1 motif (pfam00126) and a C-terminal substrate-binding 
domain (PBP2_LTTR_aromatics_like; cd08414). The gene 
orfl 7162 was renamed tciR (transfer competence inducer regulator) 
in anticipation of the results described further below. Reverse 
transcription of P. putida UWCl (ICEc/c) RNA isolated from 
exponential phase-grown cells, followed by specific PCR ampli- 
fication confirmed that the three genes are transcribed on the same 
mRNA, which ends downstream of tciR (Figure S2). This implies 
that mfsR-marR-tciR form a single polycistronic unit. 

Effects of tciR, marR, and mfsR deletions on ICEc/c transfer 

In order to more precisely investigate the role of the three 
regulators on ICEc/c transfer, their open reading frames were each 
individually and partially deleted in separate strains, namely P. 
putida UWCl (ICEc/c- Am/iP, strain 4322), UWCl (ICEdc-AmarR, 
strain 4372), UWCl {IGEdc-AmfsR-A'marR, strain 3453) and P. 
putida UWCl (ICEc/c-AtoP, strain 4321) (Figure IB, Table 1). 




2738 2961 2962 3543 
ICEclc.WT -I- - - - 
mfsR::Km - 4- -i- - 
^mfsR-'marR- - - + 




2737 4322 4646 4321 4649 4372 4804 



ICEclc.WT -I- ----- - 

Amfe/? - -I- -I- - - - 

MciR --- + + -- 

AmarR - - - - . + + 

Tn7-mfe/? - - + - - - 

Tn7-ta7? ---- + -- 
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Figure 2. ICEc/c transfer frequencies from P. putida UWC1 
donors with different ICEc/c genotypes. (A) and (B), Independently 
carried out transfer experiments using the Indicated strain sets. Bars 
show mean transfer frequencies as transconjugant colony forming units 
(CFU/ml, growing on 3CBA, Km- or Gm-resistant) per donor CFU/ml 
from biological triplicates, and the corresponding standard deviations. 
Letters above bars Indicate statistically significantly different groups per 
panel In an Anova with post hoc Tukey-Kramer test (P<0.001), with the 
same letter pointing to the absence of statistically significant 
differences. 

doi:1 0.1 371/journal.pgen.1 004441 .g002 

ICEc/c transfer frequencies in plate-mating assays with a gentami- 
cin-resistant P. putida UWC 1 as recipient were 2 • 1 0^-fold lower for 
UWCl donors with ICEc/c having an internal deletion in tdR 
compared to intact ICEc/c (Figure 2B). Complementation of the 
ICEdc- AtdR mutation with a single copy mini-Tn7 transposed 
fragment containing the tdR gene under the Pn^sR'Promoter 
(strain 4649, Figure IC) restored transfer, even to much higher 
levels than wild-type ICEc/c (Figure 2B). 

ICEc/c transfer frequencies were 27-fold lower for UWCl 
donors with ICEc/c having an internal deletion in marR compared 
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Table 1. Strains used in this study and tlieir specifications. 





strain name 


Strsin collection 
number 


Relevant characteristics 


Reference or 
source 


Escherichia coti DH5a 






[36] 


£ coli DH5a kpir 






V. de Lorenzo 


E. coli BW20767/pRL27 


1853 


tra^, pRL27 containing hyperactive mini-Tn^ element (or/V, Km"). 


[31] 


Pseudomonas !<nacl<mussii B13 


78 


Original host of ICEc/c {2 identical copies). 


[43] 


Pseudomonas putida UWC1 


1291 


plasmid-free derivative of P. putida KT2440, Rif*^ 


[44] 


P putida UWC1 (Nal) 




Spontaneous Nal'^-mutant of 1291. 


This study 


P putida UWCGC 


2744 


Single copy mini-TnZ-Ptac-echer/y insertion, Gm'^ 


[33] 


P putida UWCl 


2756 


Single copy mini-Tn5-jimX-gfp insertion, Km"^ 


[45] 


P. putida UWCl (ICEdc) 


2737 


Derivative of strain 1291 with one ICEc/c copy integrated into tRNA^'^-5. 


[45] 


p. putida UWCl (ICEc/c) 


2738 


As 2737, but integrated into tRNA^'^-6. 


[45] 


P. put/do UWCl {ICEc;c-Km",9033) 


2961 


Transposon mutant of strain 2737 with a Km'^-gene inserted at 
nucleotide position 19033 in ICEc/c. 


This study 


P putida UWCl (ICEdc-Km",86,8) 


2962 


Transposon mutant of strain 2737 with a Km'' -gene inserted at 
nucleotide position 18618 in ICEc/c. 


This study 


P put/do UWCl (ICEdc-Amftfi-A'marfi) 


3453 


Derivative of strain 2737 with mfsR and part of marR deleted (from 
nucleotide position 18395 to 19166). 


This study 


P puf/do UWCl (ICEdc) 
+ Pmt-gfp/P,nR-echerry 


3531, 3532, 3533 


Derivatives of strain 2737 with single copy random insertion of a 
mini-Tn-Pint-g/p/PinR"^'^^^''''/' ^'^^ 


This study 


P. putida UWCl (ICEdc-Afdfi) 


4321 


Derivative of strain 2737 with an internal deletion in tclR 
(from nucleotide position 17164 to 17985). 


This study 


P putida UWCl (ICEdc-Amftfi) 


4322 


Derivative of strain 2737 with an internal deletion in mfsR (from 
nucleotide position 18581 to 19143). 


This study 


r. putida UWCl (ICEc/c-A mam) 


4372 


Derivative of strain 2737 with an internal deletion in marR (from 
nucleotide position 18032 to 18468). 


This study 


R putida UWCl (ICEc/c-AmfsR) 
+ Pint-g^p/PinR-ec/ierry 


4469, 4470, 4471 


Derivatives of strain 4322 with single copy random insertion of a 
mini-Tn-Pint-g^p/PinR"^'^^^''''/' Km"* 


This study 


P. putida UWCl (ICEc/c-Amo/"/?} 
+ Pint-g^p/PinR-ec/ierry 


4475, 4476, 4477 


Derivatives of strain 4372 with single copy random insertion of a 
mini-Tn-Pint-g/p/PinR"^'^^^^''/' Km'' 


This study 


p. putida UWCl (ICEdc-Afc/K) 
+ Pint-gfp/PmR-echerry 


4479, 4480, 4481 


Derivatives of strain 4321 with single copy random insertion of a 
mini-Tn-Pjnt-g/p/PinR"^'^^^''''/' Km"* 


This study 


P. putida UWCl (ICEdc-Amftfl, gfp) 


4612 


Derivative of 4322 having a gfp gene inserted downstream of intB13. 


This study 


P. putida UWCl (ICEc/c-AmrsK, 
mini-Tn7:mft/?) 


4646 


Derivative of 4322 carrying a single copy mini-Tn7 insertion of 
the mfsR gene under its own promoter. 


This study 


F. putida UWCl (ICEc/c-Afc/n, 
m i n i-Tn 7:mfsR^'-marR^'-tciR) 


4649 


Derivative of 4321 carrying a single copy mini-Tn7 insertion of the (frameshifted) 
mfsR and marR genes, plus the intact tclR gene under the mfsR promoter. 


This study 


r. putida UWCl (ICEc/c-Amom, 
m i n i-Tn 7:mfsR^'-marR-tciR^') 


4804 


Derivative of 4372 carrying a single copy mini-Tn7 insertion of the (frameshifted) 
mfsR and tciR genes, plus the intact marR gene under the mfsR promoter. 


This study 


0 pufida UWCl rnini-TnS-Pp-ifjR-mc/if/'rj/ 


3482 


Single copy insertion of a mini-Tn5 mfsR promoter-mc/ie/ry fusion, Km'^ 


This study 


P. putida UWCl mini-Tn7-mfsR 
mini-Tn5-PmfsR-^c^e''ry 


4302 


Derivative of strain 3482 but with a mini-Tn7 insertion containing the 
intact mfsR gene expressed from its own promoter. Km'', Gm'^ 


This study 


P. putida UWCl (ICEc/c), 
mini-Tn5-PmfsR-^'^^^''ry 


3497 


Derivative of 2737, single copy insertion of a mini-Tn5 mfsR promoter-mc/ierry 
fusion, Km"* 


This study 


P. putida UWCl (ICEc/c-Amfsff), 
m i n i-Tn5-PmfsR-mc/ie/Ty 


3606 


Derivative of 3453, single copy insertion of a mini-Tn5 mfsR promoter-mc/ierry 
fusion. Km'' 


This study 


P. puf/do UWCl (\CEclc-AmfsR), 
mini-TnZ-mft/?, m\n'i-Tr\5-Pmf sR-mcherry 


4282 


Derivative of 3606, but with a mini-Tn7 insertion containing the intact mfsR 
gene expressed from its own promoter, Km'', Gm" 


This study 



doi:l 0.1 371 /journal.pgen.l 004441 .tOOl 



to intact ICEc/c (Figure 2B). Complementation of the ICEc/c- 
AmarR mutation with a similar single copy mini-Tn /-maii? 
insertion did not change transfer rates (strain 4804, Figure 2B). 
This suggests that the effect of the marR deletion on ICEc/c transfer 
is rather due to polar disturbance of the expression of the 
downstream-located tciR. 

In contrast, ICEc/c elements with mfsR deletions [i.e., P. 
putida UWCl {ICEdc-AmfsR) and UWCl {ICEclc-AmfsR-A'marR}] 



transferred with 25- and 15-fold higher frequencies than wild type 
ICEc/c, respectively (Figure 2A, B). Complementation of the 
ICEclc-AmfsR mutation with a single copy mini-transposed mfsR 
gene under control of its own promoter reduced ICEc/c-Am/si? 
transfer frequencies by 10*-fold, also here much stronger than 
predicted from wild-type ICEc/c itself (Figure 2B). These results 
suggested that tciR is the actual regulator of ICEc/c transfer, and 
further that mfsR is regulating expression of the mfsR-marR-tciR 
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operon. Since MfsR is expected to be a repressor, its deletion 
would lead to higher expression of the downstream genes marR 
and kiR, which results in increased ICEcfc transfer. The effect of 
the transposon insertions in mfsR (i.e., lower ICEcfc transfer rates. 
Figure 2A) seems therefore due to a polar effect on marR-tciR 
expression. 

tdR encodes a global activator of the genes in the ICEc/c 
core region 

Next, we examined the effect of regulator)- gene deletions on 
gene expression of ICEc/c as a whole, using semi-tiling microarray 
analyses (Figure 3). When P. putida UWCl (ICEcfc) wild-type cells 
are growing exponentially on 3-chlorobenzoate (3CBA), expres- 
sion from the genes in the ICEc/c core region (roughly the second 
half of ICEf/f) plus the integrase intBl'S is silent, whereas they are 
highly transcribed when cells are in stationary phase (Figure 3A). 
Among others, the core region encodes genes implicated in ICEc/c 
conjugative transfer [32,33]. P. putida with mutant ICEcfc lacking 
either marR or tciR strongly diminished expression in the core 
region and of the integrase gene in stationary phase when 
compared to wild type (Figure 3C, Figure S3). Lower core and 
integrase gene expression explains the lower ICEc/c transfer rates 



from these mutants (Figure 2B). In contrast, mfsR deletion resulted 
in much higher expression from the ICEcfc core genes in 
exponentially growing cells (Figure 3B), and even slighdy higher 
expression in stationary phase than in wild-type ICEc/c (Figure S4), 
which explains the 10- to 100-fold ICEt/c higher transfer rates 
(Figure 2B). Expression of the mfsR-marR-kiR cluster itself was the 
same in the tciR and marR deletion mutants, and no different to the 
wild-type (Figure 4A, C, E). In contrast, expression of the rnfsR- 
marR-tciR cluster was higher in the mfsR deletion mutants than in 
wild-type, both in exponential and stationary phase cells 
(Figure 4B, D). Since gene expression from ICEc/c is similar in 
mutants lacking rnfsR alone or mfsR plus the first 117 bp of marR 
(Figure 84), we conclude that it is the LysR-type regulator encoded 
by teiR, which is the main activator for ICEcfc core gene 
expression. 

Microarray analysis also helped to understand the behaviour of 
the mfsR Km-insertion mutant (Figure 3D). As for the kiR deletion 
mutant, expression of the ICEcfc core region and of the integrase 
was dramatically lower than wild-type in stationary phase cells 
(Figure 3D, Figure S4). On the other hand, both mfsR deletion and 
mfsR Km-insertion mutants showed intxeased expression of a 
group of genes on IGEe/c coding for a putative efflux system 
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Figure 3. Differential expression of the ICEc/c gene region from micro-array data in selected mutant ICEc/c versus wild type in P. 
putida UWCl. (A) Differential expression of the ICEc/c region between stationary and exponential phase cells of wild type P. putida UWCl (ICEc/c). (B) 
Differential expression of the ICEc/c region between the mfsR deletion mutant and wild-type, in exponentially growing cells. (C) Comparison of the 
tciR deletion mutant and wild-type, in stationary phase cells. (D) Comparison of the mfsR transposon insertion mutant versus wild-type, in stationary 
phase cells. Dots indicate the ^log-fold change of hybridization signal per microarray probe in the comparison, plotted at their distance along the 
ICEc/c sequence (X-axis; in kb). Regions of interest on ICEc/c are redrawn as grey boxes at the bottom of each section (+ or - indicate the DNA strand 
on which the region is encoded). Separate displays indicate expression differences on the plus- (open symbols) or the minus-strand (closed symbols). 
Grey bars in the background indicate the two-fold cut-off level. For a complete set of microarray results, see Figures S3 and S4. 
doi:1 0.1 371/journal.pgen.1 004441 .gOOB 
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Figure 4. Detailed view on the differential expression of the mfsR operon in P. putida ICEc/c wild-type or mutants. (A) MarR deletion 
mutant versus wild-type. (B) mfsR deletion mutant versus wild type. (C) tciR deletion mutant versus wild-type. (D) mfsR-'marR deletion mutant versus 
wild-type. (E) marR versus tciR deletions. (F) mfsR-transposon insertion mutant versus wild-type. Panels show ^log-fold change of expression level per 
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dol:1 0.1 371/journal.pgen.1 004441 .g004 



(Figure 3B, D). Detailed inspection oi mfsR-marR-tciR expression in 
the Km-insertion mutant revealed that the first 160 bp of mfsR, 
upstream of the Km'^-gene insertion were higher expressed than in 
wild-type cells (Figure 4F). In contrast, the downstream genes marR 
and tciR were lower expressed compared to wild-type and to the 
mfsR deletion mutant (Figure 4B, D, F). This confirmed, therefore, 
that insertion of the Km^-gene had caused a polar effect on 
expression of marR and tciR, which explains the strongly 
diminished expression of the ICEcfc core genes in stationary phase 
in the mfsR Km-insertion mutant, and decreased ICEc/c transfer. 

Inserting the presumed mfsR promoter region upstream of a 
promoterless mcheny gene in single copy on the chromosome of P. 



putida UWCl without ICEc/c produced strong and homogenous 
mCherry expression among all cells (Figure 5A, strain 3482). 
Inserting into this strain a single copy mfsR gene expressed from its 
own promoter abolished mCherr)' expression (Figure 5A, strain 
4302). Expression of mCherry from PmfsR in P- putida UWCl 
(ICEc/c) was very low, whereas disruption of mfsR on ICEc/c again 
resulted in high mCherry expression (Figure 5A, strain 3606). 
Complementation of this strain by a single copy mfsR gene under 
its own promoter caused repression of mCherry expression 
(Figure 5A). All these data are consistent with the hypothesis that 
MfsR is repressing expression of itself and the downstream located 
marR and tofl genes. 
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Figure 5. Effect of mutations in tlie mfsR region on tlie expression of the PmfsR-, Pint- and Pi„R-promoters of ICEc/rin P. put/da UWC1. 

(A) mCherry expression from tine mfsR promoter added in single copy to the chromosome of the indicated P. putida UWC1 strains (relevant genotypes 
and strain numbers specified below the graph). mCherry expression measured on individual cells (n = 1000) by epifluorescence microscopy in late 
exponential phase of cultures grown on 10 mM succinate and expressed as box plots (AU, arbitrary units at 20 ms exposure time). (B) Scatter plot of 
GFP and mcherry fluorescence in single cells of P. putida UWC1 (ICEc/c) wild-type, AmfsR, AtciR or AmarR deletions, equipped with a single copy mini- 
transposon containing the Pint-egfp and PjnR-mc/ierjy fusions. Panels show expression of both markers at different growth phases, as indicated, with 
colors representing genotypes with independent mini-Tn5 insertions. Note as example the subpopulation of wild-type cells (dotted ellips) expressing 
both reporters, compared to the majority of cells in the mfsR deletion mutant but a complete absence of such subpopulation in the tciR and marR 
deletion mutants. 

doi:1 0.1 371 /journal.pgen.1 004441. g005 



mfsR deletion leads to an increase in the number of cells 
activating ICEc/c 

We then tested whether changed ICEc/c transfer rates and core 
gene expression were in fact due to changes in the proportion of 
cells activating ICEc/c. Hereto, a double promoter-reporter 
construct, carrying FintSfP ^"^d ¥^^-echeny was inserted in single 
copy on the chromosomes of P. putida UWC 1 (ICEc/c), P. putida 
UWCl {ICEclc-AmfiR), P. putida UWCl (ICEc/c-A'marii) and P. 
putida UWCl {ICEclc-AtciR) (Table 1). Pin, and Pi„R are the 
respective promoters for the integrase gene intBlS (integrated 
form) and the integrase activator gene inrR. Previous studies 
showed that both promoters are active only in a small 
subpopulation of cells, which are representative for transfer 
competent cells [26,29,30]. Consistent with previous data, the 
subpopulation of P. putida UWCl (ICEc/c) wild-type cells 
expressing P;„, and PjnR in stationary phase suspended cultures 
represented a few percent (Figure 5B). In contrast, deletion in tnfiR 
resulted in 80-100% of cells expressing Pint- and PinR-promoters 



(Figure 5B). Expression of both promoters in P. putida UWCl 
{ICEclc-AnifsR) occurred in early stationary phase whereas in wild- 
type cells their expression is maximal in late stationary phase 
(Figwe 5B). Conversely, P. putida UWCl {IGEclc- AtciR) and P. 
putida UWC 1 (ICEctc-AmarR) did not produce any detectable Pinr 
or PinR-expressing cells, neither in exponential nor in stationary 
phase (Figure 5B). Considering a detection limit by microscopy of 
~1 fluorescent cell among 1000-10,000 non fluorescent cells, the 
absence of detectable Pint- or PinR-expressing cells in those mutants 
would be in accordance with absence of ICEc/c core gene 
activation on microarrays (Figure 3D) and lower transfer 
frequencies (Figure 2B). 

Mutants with mfsR deletion in ICEc/c face a strong fitness 
cost 

Given that P. putida UWC 1 carrying ICEclc-ArrifsR transferred at 
a much higher rate than wild-type ICEc/c, and also expressed both 
Pint- and PinR- promoters in almost aU cells, we wondered why 
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Table 2. Effects of the mfsR deletion on the growth characteristics of P. putida UWCl carrying ICEc/c. 





strain 




Generation time (min)^ 


Survival rate (%)^ 




Name 


Number 


MM Succ 


MM 3CBA 


MM Succ 


MM 3CBA 


P. putida UWCl (ICEdc) 


2737 


74.4 ±6.8^ 


186±11 


100±12 


113±8.5 


P. putida UWCl (ICEc/c-Amftft) 


4322 


68.2±8.1 


184±3 


0.3 ±0.3 


5.6±2.5 






p = 0.38'' 


p = 0.77 


p = 0.0001 2 


p = 2.96- 10"^ 



^'Generation time was calculated as In2/[i, whereby \i (min'^) is the slope of the regression line on a plot of the loQcuiture turbidity versus time from at least 5 points during 
exponential phase. Coefficients of determination (R^) were >0.96 for each growth curve. 

^'Survival rate in stationary phase of cultures on the indicated media was calculated as the ratio of the number of CFU/ml counted on MM+3CBA and the number of 
CFU/ml on MM+succinate plates. Succ, succinate; 3CBA, 3-chlorobenzoate. 
^'Calculated standard deviation from triplicate measurements. 
'"Calculated p-value in a two-tailed Student's t-Test using equal variance. 
doi:l 0.1 371 /journal.pgen.l 004441 .1002 



such mutants did not become selected spontaneously. Both P. 
putida UWCl (IGEc/c) wild-type and (ICEcfc-Am/yi?) displayed 
statistically indistinguishable generation times during exponential 
growth on minimal medium with either 3CBA or succinate as 
carbon source (Table 2), although P. putida UWC I (ICEc/c -Am/si?) 
went through a longer lag phase (Figure S5A). In contrast, the 
proportion of colony forming units (CFU) in samples taken from 
stationary phase cultures on 3CBA or succinate and plated on 
3CBA solid medium was dramatically reduced for P. putida UWC 1 
(ICRclc-tSmfsR) (Figure 6). Whereas the number of colonies formed 
on 3CBA plates was the same as that on succinate for P. putida 
UWCl (ICEc/c), only 0.3-5.6% appeared on 3CBA plates for 
UWCl {lCEclc-^mfsR) (Table 2). Also in absolute terms, the 
number of CFU/ml for UWCl (ICEc/c-Am/ri?) ceUs taken from 
stationary phase cultures both on 3CBA and succinate was lower 
than that for UWCl (ICEc/c) (Figxire 6). Moreover, 8 of 10 tested 



colonies of UWCl (ICEck-AmfsR) grown on MM plates with 
succinate did no longer amplify the clcA gene of ICEc/c (not 
shown), the remaining two still being able to grow on 3CBA. 
Furthermore, half or more of microcolonies formed from UWC 1 
[ICEclc-AmfsR) with a single copy Fint-fgfP insertion showed 
incidence of malformations and cell lysis, similar to what was 
reported previously for nutrient-reactivated tc cells [30] 
(Figure 6B), but cells in the other microcolonies divided with 
generation times even slightly faster (1.49±0.15 h) than those in 
microcolonies of P. putida UWCl (ICEc/c, 1.79±0.09 h). This 
indicates that the mfsR deletion in ICEc/c imposes a strong fitness 
cost on P. putida UWC 1 . Survival of UWC 1 was restored to wild- 
type level when the ICEc/c-Am/si? was complemented by the mini- 
Tn7 inserted mfsR gene (Figure S6). In contrast, neither P. putida 
UWCl (ICEc/c) with marR or tciR deletion, nor the mini-Tn7 
complemented strains of P. putida UWCl (ICEc/c-A/c/i?) and 



A B 




Figure 6. Fitness loss of P. putida UWCl (ICEc/c) caused by the mfe/? deletion. (A) Survival of P. putida UWCl (ICEc/c, 2737) and P. putida 
UWCl (ICEc/c-AmfsR, 4322) pregrown in suspended culture to stationary phase on 3-chlorobenzoate (3CBA) or succinate (succ), and plated from there 
on 3CBA or succinate agar. (E.g., 3CBA-3CBA, suspended culture on 3CBA, plated on 3CBA agar). Survival expressed as colony fornning units (CFU) on 
the agar plate per ml of stationary phase culture. Data bars indicate the average from independent biological triplicates. Error bars indicate the 
calculated standard deviation from the average. (B) Phase-contrast micrographs at 1000-fold magnification of microcolonies of P. putida UWCl 
(ICEc/c) and P. putida (ICEc/c-Am/s/? Pin,-eg/p, 4612) growing on agarose surface supplemented with 0.1 mM 3CBA. Shown are a regular stationary 
phase microcolony of P. putida with wild-type ICEc/c and a transfer competent microcolony (tcm), occurring at 1-3% frequency as reported previously 
[30]. For comparison, massive lysis (white arrows) and cellular malformations formed in many microcolonies of P. putida UWCl with the mfsH deletion. 
doi:1 0.1 371/journal.pgen.1 004441 .g006 
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Figure 7. Conservation of fc//? analogues in putative ICEc/c-like regions in a variety of other bacterial genomes. Illustration represents 
tciR analogues (identified on the basis of a BLASTN E-value lower than 1 -10"'^), and surrounding relevant gene regions in the indicated bacterial 
genomes (species name, accession numbers) compared to ICEc/c. Genes are indicated as in the respective genome accession. Rectangles show 
annotated genes and their orientation (top, orientation towards the left; bottom, gene orientation towards to right); common colors indicate similar 
predicted functions. Stippled rectangles indicate common gene regions inferred from Artemis comparison, but not present in the respective 
annotation. ICE were inferred from (i) more than 75% nucleotide identities across the complete core region of ICEc/c, and within a 1-100 kb window 
from the tciR position, and (ii) the presence of an integrase gene (in brown) within a 5-20 kb window from the fc//?-analogue. Note how some 
genomes carry multiple different ICE from the same family (e.g., Achromobacter xylosoxidans, Acidovorax sp. strain JS42), and further how pair-wise 
identical ICE regions (shaded in grey) occur between different genomes. Finally note how the fc//?-analogues often co-occur with a xer-type 
regulatory gene on the other strand (light green), and a further lysR gene member (yellow), but in none of the cases shown here with an mfsR 
counterpart (in red). 
doi:1 0.1 371/journal.pgen.1 004441 .g007 
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(ICEclc-AmarR) were impaired in survival compared to UWCl 
(ICEf/c) (Figure S6). 

tciR is a very widespread ortholog among ICE closely 
related to ICEc/c 

Using bioinformatic queries, we retrieved orthologs to tciR and 
mfsR from sequenced bacterial genomes, and examined manually 
whether they occur in chromosomal regions qualifying as ICE 
(e.g., presence of an integrase gene nearby, see Materials and 
Methods). Interestingly, orthologs to the individual components of 
the mfsR operon are widespread but rarely occur in the same 
configuration (Figure 7). So far, the ICEc/c mfsR-marR-tciR 
configuration is only found in the ICEc/c variant of Burkholderia 
xenovomns LB400, whereas the Tn^577-element of Acidovorax sp. 
strain JS42 and Aenmonas hjdwpohila SSU (accession number 
AGWR01000022.1) both carry an mfsR homolog and a MarR- 
type regulator immediately downstream, but no tciR equivalent 
nearby. There is a tciR ortholog in Acidovorax sp. strain JS42, but 
not on the same chromosomal region as mfsR (Figure 7). 

On the other hand, tciR seems much more widespread among 
ICEs, as homologs can be found in ICEc/c-like elements Gil and 
GI6 oi Bordetella petrii DSM 12804, in PAGI-2 of P. aeruginosa strain 
C, in diverse ICEs of X. campestris pv. vesicatoria str. 85-106', 
suspected ICEs in Herminiimonas arsenicoxydans, Cupriavidus metalli- 
durans CH34, and Tolumonas auensis DSM 9187, among several 
dozens of others (Figure 7). Given that TciR is such a common 
regulator found in ICEs of the ICEc/c family, it might be similarly 
implicated in their transfer control. This regulation is likely 
dilferent in detail from ICEc/c, given the absence of an nifsR and 
marR. 

Discussion 

ICEc/c has two distinctive modes of existence: the integrated 
form, which is transmitted vertically, and the circular form, which 
can be horizontally transferred. Previous work in our laboratory 
has shown that the transition between these two states occurs in 
only a few percent of cells in a population under stationary phase 
conditions [26,30,34]. We have recentiy suggested to name cells in 
which the molecular decision occurs to activate the ICEc/c 
horizontal transfer mode transfer competent (tc) cells [30] . Single cell 
time-lapse experiments indicated that ICEc/c transfer - at least 
insofar as detectable by microscopy, only occurs from tc cells, 
which can be distinguished through simultaneous expression of 
fluorescent proteins from single copy transcriptional fusions to the 
Pjjit and PinR promoters of ICEc/c [30] . Activation of those two 
promoters is the likely outcome of a multi-step regulatory cascade 
that orchestrates expression of some fifty genes [32], but the key 
factors that determine the onset of this cascade and control the 
extent of bistabUity are still obscure. Previous work provided 
evidence for the role of the stationary phase sigma factor RpoS in 
activation of ICEc/c promoters, and we could show that tc cells on 
average have higher levels of RpoS [34] . In the present study, we 
report the discovery of a cluster of three regulatory genes, two of 
which globally control ICEc/c activation and transfer, and 
additionally maintain a cap on fitness loss induced by the ICE 
(Figure 8). 

The three regulatory genes occur in a unique configuration of a 
TetR-type repressor (encoded by mfsR), followed by a MarR-type 
{orfl7984] and a LysR-type regulator [tciR). Transcript and 
microarray analysis of the locus in wild-type and mutant ICEc/c, 
plus analysis of reporter gene expression from the mfsR promoter 
in a variety of host backgrounds, showed that the three genes are 
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Figure 8. Model for regulation of ICEc/c transfer competence. 

iVlfsR autoregulates expression of itself and of the TciR activator, 
without which ICEc/c transfer decreases by 2-10^-fold. TciR may activate 
specifically one or more promoters on ICEc/c, such as the RpoS- 
dependent PinR-promoter [29], or a promoter upstream of the parA-WWe 
gene [32], Expression of the PinR-promoter occurs preferentially in cells 
having highest RpoS levels, and only 1-3% of cells in a population in 
stationary phase visibly express reporter gene from PjnR [29]. InrR 
transmits bistable activation through an unknown process to the intBIS 
promoter [26], and possibly simultaneously to other promoters for the 
genes for the conjugative system in the ICEc/c core region [32]. The 
exact mechanism of arisal of bistability is unclear as yet. For gene 
locations on ICEc/c, see Figure 1A. 
doi:1 0.1 371/journal.pgen.1 004441. g008 



expressed as a polycistronic unit and are under autoregulatory 
control by MfsR (Figure 5A). 

Precise gene deletions and complementations indicated that tciR 
is likely the main global regulator of ICEc/c transfer activation. 
Deletion of tciR caused a 2-10^-fold lower frequency of ICEc/c 
transfer compared to wild type (Figure 2B), silenced expression of 
the ICEc/c core region (Figure 3B) and reduced the proportion of 
cells expressing Pju, and Pi„R in stationary phase (Figure 5B). 
Complementation of the tciR deletion on ICEc/c with a single copy 
mini-Tn7 inserted mfsR!'' -marR!' -tciR fragment (Figure IC) fused to 
the mfsR promoter region restored the expected phenotype 
(Figure 2B). TciR may act either directly as a regulator on a 
variety of individual ICEc/c core promoters, or as a "master" 
regulator in a hierarchical activation cascade (Figure 8). 

The role of marR is less clear and as yet unsolved. Deletion in 
marR resulted in essentially the same ICEc/c transcriptome profile 
as deletion in tciR (Figure 3D). It also resulted in a lower transfer 
frequency than wild-type but not as low as the deletion in tciR 
(Figure 2B), and produced no detectable reporter gene expression 
from P;,j, or Pj^^ (Figure 5). In contrast, complementation of the 
marR deletion on ICEc/c by a single copy marR gene through mini- 
Tn7 delivery (Figure IC) did not restore ICEc/c transfer 
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(Figure 2B). Furthermore, mutants with deletions in mfsR or mfsR 
plus the first 1 16 bp oi marR (Figure IB), behaved quite similar in 
transfer fi-equency (Figure 2) and showed similar ICEcfc tran- 
scriptomes (Figure S4). The marR gene therefore seems to have no 
clear role in ICEc/c core gene expression. 

The most surprising effect of deletions in mfsR was a complete 
deregulation of the ICEcfc core gene expression. This became 
obvious from frequencies of ICEc/c transfer being 10-100 fold 
higher than wild type, approximating 1 transfer per donor cell 
(Figure 2B). The deregulation was also obvious in microarray data 
showing the ICEc/c core region in the mfsR deletion mutants being 
already transcribed in exponential phase (Figure 3B, Figure S3). 
Finally, 80-100% of individual cells in stationary phase expressed 
the reporter genes from Pij,t and PjnR in the mfsR deletion strain 
compared to 3-5% in wild-type (Figure 5B). This can be explained 
by the fact that deletion of mfsR would abolish autorepression, 
which would lead to constant high expression of marR and tciR. 
This overinitiates ICEcfc core gene expression, leads to more cells 
entering the tc state and to higher transfer rates. The balance of 
mfsR control appears to be extremely delicate, since even 
complementation with a single gene copy under control of the 
original promoter results in a stronger effect than the wild-type, 
both for tciR and mfsR (Figure 2B). The delicate balance became 
also obvious from the polar effects of insertion of the Km- 
resistance gene within mfsR, leading to decreased transcription of 
marR and tciR, diminished core gene expression (Figure 3D) and 
reduced transfer rates (Figure 2A). 

The finding that deletion and complementation of m/si? or tciR 
drastically changes the proportion of cells activating ICEcfc, could 
imply that the bistabUity seen in wild-type situation (i.e., 3-5% of 
cells in stationary phase becoming transfer competent) is a result of 
feedback at this locus. Cells activating ICEcfc in the wild-type 
situation could arise as a consequence of "sloppy" control by 
MfsR, incidentally causing a few cells to escape its control and 
transcribing rmrR and tciR. We think this is an unlikely scenario, 
because mCherry expression from the mfsR promoter is homog- 
enous among cells (Figure 5A). Alternatively, there might be a 
chemical ligand that specifically binds to MfsR in a small subset of 
cells, upon which its repression is relieved in those cells. The 
resulting TciR would then be the necessar)' activator to trigger 
ICEcfc core expression in cells with on average highest RpoS levels 
[29] (Figure 8). On the other hand, even though mfsR may be the 
first level of control, bistabUity may also originate at later 
checkpoints in the regulatory cascade, which depend on the 
presence of sufficient TciR. 

Quasi-global appearance of transfer competence across all cells 
in the mfsR deletion mutant resulted in massive fitness loss 
(Figure 6), which became evident at two levels. First of all, time- 
lapse observations indicated lysis and aberrant cell growth in more 
than 50% of microcolonies (Figure 6B). This lysis and growth 
arrest are similar to what we previously described as being a side 
conserjuence of becoming transfer competent in wild-type cells 
[30], and is caused by the parA-shi gene products on ICEcfc [30]. 
Secondly, there was a strong loss of the capacity to grow on 3CBA 
among cells sampled from stationary phase cultures of the mfsR 
mutant compared to wild-type (Figure 6A, Figure S6), indicative 
for loss of ICEcfc and counterselection against maintaining ICEcfc- 
AmfsR. However, those cells that maintained ICEcfc-Am/si? could 
still grow on 3CBA and showed indistinguishable exponential 
growth rate (Table 2). This paradox can be understood when 
modeling the number of tc cells in batch culture populations for 
ICEcfc wild-type (probability of tc arisal, Pt,., in stationary phase of 
0.025) and for the ICEclc-AmfsR mutant {P^ - 0.5). This model 
(Figure S5B) shows that whereas a large proportion of tc cells 



appear in ICEcfc-Am/ri? mutant cultures in stationary phase, these 
can only divide 2-3 times upon reinoculation into fresh medium 
before lysing. This causes an apparent prolongation of a lag phase 
visible as stagnant culture turbidity, but does not influence the 
overall predicted population exponential growth rate in batch 
culture (Figure S5B). 

As expected from the postulated role of TciR, its complemen- 
tation in trans also leads to increased ICEcfc transfer, but 
interestingly, only the mfiR deletion caused strongly decreased 
cell survival (Figure S6). We therefore hypothesize that ICEcfc 
activation may follow two separate processes: transfer and tc cell 
growth arrest [30], that may both be initiated at the rnfsR locus. 
Deleting mfsR would then deregulate both processes, whereas 
expressing tciR in trans would only increase activation through the 
transfer branch (Figure 8). 

The configuration of the rnfsR-marR-tciR operon of ICEcfc is 
unique, but tciR alone is a very common part of ICE similar to 
ICEcfc (Figure 7). We therefore speculate that mfsR-marR are a 
more recent acquisition in ICEcfc, which drastically changed the 
expression of the tciR gene. Unfortunately, expression of tfiR 
analogs in other ICEs has not been studied and very litde has been 
reported on the transferability of ICEs related to ICEcfc. The 
exceptions being GI3 of B. petrii that transfers at extremely low 
frequencies (~10'') [35], and the P. aeruginosa PAGI-2 element for 
which transfer has not been detected at all [16]. In comparison, 
wild-t}pe ICEcfc transfers at rates of 10'^ to 10''^ per donor 
(Figure 2), suggesting that it was perhaps the acquisition of the 
mfsR regulatory control that led to expression of transfer activity in 
a larger proportion of cells in the population. As we show here, the 
downside of increasing the proportion of ICEcfc tc cells is an 
increase of the proportion of cells displaying growth arrest through 
the shi-parA pathway [30]. Likely, the MfsR autoregulation evolved 
to a stage of permitting efficient transfer but avoiding too much 
fitness loss to the population. Even though the mechanistic details 
are different for ICEcfc, double control layers are more common 
for various ICEs and typically invoK-e a variety of regulators acting 
on each other and/or in response to specific chemical ligands 
[18,19,23-25]. It will be highly interesting to further study the 
mechanistic details of the control systems that maintain very low 
ICE transfer rates, and to understand whether and how such 
control can evolve to allow hyperefficient transfer. 

Materials and Methods 

Strains and culture conditions 

Table 1 lists the strains used in this study. Eschenchia coli DH5a 
(Gibco Life Technologies, Gaithersburg, Md.), E. coli DH5a Xpir, 
E. coli BW20767/pRL27 were cultured at 37°C on Luria-Bertani 
(LB) medium [36]. Pseudomonas species were cultured at 30°C on 
LB or 21C minimal medium (MM) [37] complemented with one 
of the following carbon sources: 0.5, 5, or 10 mM 3-chlorobenzo- 
ate (3CBA), 15 mM succinate or 10 mM fructose. Antibiotics were 
supplemented to the growth medium to select for maintenance of 
genetic constructions at the following concentrations: kanamycin 
(Km) 25 Hg/ml, chloramphenicol (Cm) 20 Hg/ml, rifampicin (Rif) 
50 Hg/ml, nalidixic acid (Nal) 50 |a.g/ml, gentamicin (Gm) 20 Hg/ 
ml, and ampiciUin (Ap) 100 |ig/ml. 

Strain constructions and DNA teciiniques 

DNA purification, PGR, restriction enzyme digestions, DNA 

ligations and electro-transformations were performed according to 
standard procedures [36]. Deletions in ICEcfc genes were created 
by double recombination techniques as described elsewhere 
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[29,38]. Nucleotide positions are given according to AJ617740 
(ICEc/c). Primers used for strain constructions are listed in table SI. 

For complementation of P. putida UWC 1 (ICEc&-Am/?i?, strain 
4,322) we first amplified the mfsR gene plus the 429 bp upstream 
region containing the mfsR promoter using PGR. This fragment 
was cloned into pGEM-T-easy and verified for correctness by 
DNA sequencing. The fragment containing the correct mfsR 
region was then recovered by restriction enzyme digestion with 
PstI and BamHI, and ligated into the mini-Tn7 vector pUC- 
miniTn7-Gm [39]. After transformation and verification in E. coli, 
the mini-Tn7 txjnstruct was introduced into P. putida UWGl 
(ICEclc-AmfsR) by using the pUX-BF13 helper plasmid [40]. 
Clones resistant to Gm were selected and verified by PGR for 
correct insertion of the mfsR DNA in the attTnJ locus. To 
complement P. putida UWGl (ICEcfc-Afctfl, strain 4321) and 
(ICEclc-AmarR, strain 4372) we amplified the complete mfsR-maiR- 
tciR locus including the 429-bp upstream region. This fragment 
was cloned into pGEM-T-easy and again verified for correctness 
by DNA sequencing. The fragment was recovered by digestion 
with BamHI and StuI, and ligated with the mini-Tn7 vector. A 
frameshift was then introduced in the mfiR coding region by 
digestion at the unique Ncol-site, filling in using Klenow and 
religation. This will cause premature ending of the mfsR gene 
product {mfsR^'J. A second frameshift was subsequentiy introduced 
to inactivate the marR gene product, using the unique Bsgl-site 
{marPf'). After transformation and verification in E. coli, the 
construct was introduced in P. putida UWGl (IGEcfc-Afcfi?) as 
outlined above. This procedure was repeated to create a fragment 
with frameshifts in mfsR and in tciR (using the unique Kpnl site), 
but maintaining an intact marR. This construct was introduced 
into P. putida UWGl (ICEck-AmarR). Gm-resistant clones were 
verified by PGR for the correct insertion at the attTnl-ske, and for 
the presence of IGEc/r . 

A 656-bp region upstream of mfiR was amplified by PGR and 
fused to a promoterless mcherry gene. This fragment was introduced 
in single copy on the chromosome of P. putida UWG 1 , P. putida 
UWGl (IGEc/c) or P. putida UWGl (IGEcfc-Am/ji?) using mini-TnJ 
delivery. Three independent Km-resistant colonies were verified 
by PGR for the correct insertion and stored individually. P. putida 
UWGl mini-Tn,5-P^ftR-mcAm7 and P. putida UWGl (IGEcfc- 
AmfsR) mini-Tn^-PnjfsjR-meferrc were then further used as recipient 
to introduce the mim-Tn7-mfsR construct. 

Random mutagenesis and screening 

Random mini-transposon insertions in P. knackmussii B13 were 
generated by mobilization of the pRL2 7 suicide plasmid from E. 
coli BW20767 in a biparental mating. Hereto both strains were 
each cultured overnight in 3 ml LB, pelleted down, resuspended in 
50 |xl sterile saline solution (0.9% NaGl), mixed in a 1:1 [ii/v) ratio 
and incubated on the surface of an LB agar plate for 24 hours at 
30°G. The mixture was then resuspended with 1 ml saline 
solution, which was inoculated in 100 ml MM with 0.5 mM 
3GBA plus Km to select for the mini-transposon insertion and Gm 
to counterselect against E. coli, and incubated at 30°G for 16 h 
with orbital shaking (180 rpm). An aliquot of 3 ml of this pool of 
enriched Km^ B13 mutants was used en masse as donor in a 
subsequent mating procedure. Hereto, cells from the 3 ml 
suspension were pelleted by centrifugation, washed with 3 ml 
sterile saline and mixed with 3 ml of suspension of P. putida UWG 1 
recipient, that had been grown for 16 h on LB, was pelleted by 
centrifugation and resuspended in sterile saline. The mating 
mixture was again centrifuged, the cell pellet was resuspended in 
50 |J,1 sterile saline solution and spotted on the surface of a MM 
agar plate containing 0.5 mM 3GBA. The mixture was incubated 



for 72 hours at 30°G, after which the cells were washed from the 
plate with 1 ml sterile saline, which was further serially diluted and 
plated on MM agar plates with 5 mM 3GBA plus Km and Rif to 
select for transconjugants carrying mutant IGEc/c. Individual 
colonies were purified, recultured in organized 96-well format and 
stored at -80°G after addition of and mixing with glycerol to 15% 
(»/»). Libraries were repUcated and regrown in 100 jti LB plus Rif 
for 16 h in 96-weIl microtiter plates, mixed with 100 nl P. putida 
UWGl Nal^ recipient suspension, and incubated at 30°G for 
48 h. Then 50 ixl of each well was reinoculated into 170 |a,l of MM 
containing 5 mM 3GBA plus Km, Rif and Nal, and growth was 
measured by continuous OD-measurements in a multiplate reader 
(FluoStar Omega, BMG labtech). Absence of growth was taken as 
indication for absence of IGEc/c transfer, in which case the donor 
culture was recovered for mapping of the transposon insertion. 

Insertion mappings 

DIG-labeled primers 070934 or 070935, annealing to one of the 
ends of the Km^ insert but facing outward, were used (separately) 
in single-primer PGR with DNA from mutant UWGl donors as 
templates. The reactions produced oligonucleotide probes with the 
5 '-DIG label, the serjuence of the end of the Km'*^ gene and the 
adjacent sc(]ueii('e of the IGEf/c insert position. Such products 
were used for rough localization of the insertion position by 
hybridizing to macroblot membranes (Eurogentec, UK), whose set 
of 55-mer oligonucleotides covers most of the IGEcfc genes. 
Hybridization and detection of the DIG-marker were carried out 
according to the manufacturer's instructions (Roche Diagnostics 
GmbH, Mannheim, Germany). Once the insertion was roughly 
mapped on IGEefc, PCR-based sequencing was used to determine 
exact position of the Km'^-gene insertion. 

ICEc/c transfer assays 

The frerjuency of IGEc/c transfer was determined in experi- 
mental conditions described previously [33]. P. putida UWGl 
IGEc/c wild-type or mutant derivatives were used as donors, 
whereas P. putida UWGGG (constitutively fluorescent, Gm^) or P. 
putida UWGl Km'^ were used as the recipient (Table 1). Briefly, 
donors and recipient were each cultivated on 5 mM 3GBA MM 
and 10 mM fructose MM, respectively, and combined on 0.5 mM 
3GBA agar plates as a single concentrated peUet. After 48 hours 
incubation at 30°G, mating mixes were resuspended, diluted and 
plated on 5 mM 3GBA MM agar (counting of donor GFU) or 
5 mM 3GBA Gm or Km agar (counting of transconjugant GFU). 
Transconjugants were checked by PGR and frequencies were 
expressed as the number of transconjugant GFU per donor GFU. 
Donor survival was used for the data shown in Figure S6. 

ICEc/c transcriptome analysis by microarrays 

The IGEc/c transcriptomes of P. putida UWG 1 (IGEc/c), P. putida 
UWGl (ICEc/c), P. putida UWGl (ICEc/c-Km'^igoss), P. putida 
UWGl (ICEc/c-A'mflrP), P. putida UWGl [IGEclc-A'tciR), P. putida 
UWGl {ICEclc-AmfsR-A'marR), P. putida UWGl {ICEclc-AmfsR) 
were investigated by microarray analysis, as described previously 
[32] . Total RNA was extracted from ceUs grown on 1 0 mM 3GBA 
MM, and harvested at mid exponential phase (OD(,()() = 0.6) and 
48 h after entrance in stationary phase. Reverse trans( ription 
using cyanine-dCTP among the dNTPs produced labeled cDNA 
that was further purified and hybridized on 8x15 K microarray 
slides (Agilent, Santa Glara, GA, USA). Slides were washed and 
scanned according to manufacturer's instructions (Agilent). Data 
were recovered and analyzed using GeneSpring GX. Microarray 
data can be accessed from the GEO database (accession number: 
GSE51391). 
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Time-lapse microscopy 

P. putida UWG 1 strains were precultured for 16 h at 30°C in LB 
medium, after which 100 ^ll were transferred to 20 ml fresh MM 
4 mM 3CBA medium in presence of the appropriate antibiotics. 
This culture was incubated for 96 hours at 30°C and 200 rpm 
shaking, after which the cells were 1 00-fold diluted in MM without 
C-source and inoculated on agarose surfaces {gel patches) for time- 
lapse microscopy [41]. Medium for gel patches consisted of 1% 
agarose dissolved by heating into MM with 0. 1 mM 3GBA. Gel 
patches were created by pipetting 130 |J.l of the agarose-MM- 
3GBA solution kept at ,55 °G on the surface of a circular cover 
glass (42 mm o and 0.17-mm thick), placed in an autoclaved 
perfusion chamber (POG, H. Saur, Reutlingen, Germany), 
separated with a 0.5 mm thick sUicon spacer ring and covering 
them with a second cover glass. After solidification of the agarose, 
the upper cover slip was removed and 6 (tl of the diluted cell 
suspension was placed onto the agarose gel patch. As soon as the 
drops were dried on the surface, the patches were turnc'd upside 
down and placed bacteria-facing-down on a new round cover glass 
[41]. A second silicon spacer ring was added to allow air 
circulation within the closed chamber and the glass sandwich was 
frxed into the metal cast POG chamber with a metal ring. Up to 
four patches could be placed simultaneously within a single glass 
sandwich in a POG chamber. 

Microcolony development was followed directly on a Nikon 
Inverted Microscope Eclipse Ti-E, equipped with a Perfect Focus 
System (PES), pE-100 GoolLED and a Plan Apo X, 100x1.45 Oil 
objective (Nikon), installed in a controlled temperature room 
(22°G). Ten random regions of every patch were imaged 
automatically during 48 hours with intervals of 1 h, in Phase 
Gontrast mode (10 ms exposure), eGFP (500 ms) and cGherry 
(500 ms). Images were recorded using Micro-Manager 1.4 (http:// 
www.micro-manager.org/) and fluorescence values were extracted 
using MetaMorph (Series 7.5, MDS, Analytical Technologies). 

P. putida UWC1 ICEc/c and mutant fitness tests 

Triplicates of strains UWGl gGEcZc) and UWGl IJGEclc-^mfsR) 
were grown for 16 h in LB medium at 30°G. Both strains were 
then 500-fold diluted (starting ODroo 0.001) in MM with 5 mM 
3GBA or 10 mM succinate. Upon reaching early stationary phase, 
strains were again diluted into fresh MM (starting OD 0.001) with 
the same carbon source, and growth was followed by frequent 
culture turbidity measurements (ODgoo). 24 h after reaching 
stationary phase, each replicate culture was serially diluted in 
MM and plated onto MM agar plates with 5 mM 3GBA or with 
10 mM succinate. The number of CFU/ml was scored and the 
ratio was calculated b(;twcx'n the number of GFU/ml on MM agar 
with 5 mM 3GBA and the number of GFU/ml on MM with 
10 mM succinate. 

Ten randomly chosen colonies of UWGl (IGEcfc-Am/ri?) 
cultivated in MM with succinate and grown on MM-succinate 
agar plates wx-rc rctc'sted for growth on MM agar with 5 mM 
3GBA. The presence of IGEc/c was determined by colony PGR on 
the same colonies by amplifying the dcA gene, which is carried by 
IGEcZc and the gene product of which is essential for 3GBA 
metabolism. 

Bioinformatic screening for ICE related to ICEc/c 

Homologues to tciR of IGEc/c were detected by BLASTN to the 
nr/nt database at E-value <1-10''*. The corresponding whole or 
draft genome sequences were retrieved and compaired by aligning 
to IGEc/c (Accession number AJ617740.2) using Megablast. 
Detected regions were manually recovered and searched for the 
kiR homologue and an intBlS homologue within a 1-100 kb 



window. If annotated, the presence of a gene for tRNA-Gly nearby 
the intB13 homologue was scored. Regions covering all criteria 
(i.e., homology to IGEcfc core region, kiR homologue and presence 
of integrase gene) were retained as containing putative IGE. 
Selected regions were further individually pair-wise compared by 
using the Artemis Gomparison Tool within the WebAGT service 
[42]. 

Supporting Information 

Figure SI Outline of the random mutagenesis and subsequent 
selection procedure. Original IGEc/c-host P. kmckmussii B13 is 
randomly mutagenized by miniTn5-mediated insertions of 
kanamycin resistance inserts (Km^). Mutant BIS are selected by 
culturing cells on minimal medium (MM) with Km and 3- 
chlorocatechol (3GBA) as sole carbon and energy source. The pool 
of B 1 3 mutants is cultured in batch and mixed with recipient strain 
P. putida UWGl (resistant to rifampicin, Rii*). The mixture is 
incubated in mating conditions for 72 hours and plated on MM 
with 3GBA, Km and Rif, to select for transconjugants. Individual 
colonies of transconjugants were restriked and organized into a 
mutant library in 96-wen plates. Each mutant is used as donor in a 
new 96-well mating with recipient P. putida UWG 1 resistant to Rif 
and nalidixic acid (Nal'^), Individual mating mixtures are grown 
on MM agar with 3GBA, Km, Rif, Nal in order to select for 
transconjugants. In absence of transconjugant growth, the donor 
of that particular mating was recovered from the library and had 
its insertion position mapped. 
(TIF) 

Figure S2 Reverse transcriptase polymerase chain reaction 
analysis of transcription in the mfsR-marR-tciR region. (A) 
Amplification of specific regions on reverse-transcribed (+) or not 
(-) mRNA purified from exponentially growing P. putida UWGl 
(IGEc/c) cultures on 3GBA, compared to amplification on purified 
DNA. (B) Schematic overview of the location of the used primers 
for the reverse transcription reaction (RT) and for the amplifica- 
tion of the gene regions. 
(TIF) 

Figure S3 Pair-wise comparisons of expression in the ICEck 
area by microarray analysis. (A) P. putida UWGl (IClLck-AmarR, 
strain 4372) compared to P. putida UWG 1 (ICEc/c, strain 2737). (B) 
P. putida UWGl {ICEclc-AtaR, strain 4321) versus P. putida UWGl 
(ICEcfc, strain 2737), (C) P. putida UWGl (ICEctc-AmarR, strain 
4372) versus P. putida UWGl {ICEclc-AtciR, strain 4321), Panels 
indicate comparisons of exponentially growing or stationary phase 
cells, with hybridization signals on the plus- (open symbols) or 
minus-strand (closed symbols) of IGEc/c, Dots indicate the ^log-fold 
change of hybridization signal per microarray probe in the 
comparison, plotted at their distance along the IGEc/c sequence 
(X-axis; in kb). A scheme of IGEcZc is redrawn at the bottom of 
each section, with regions of interest as grey boxes (+ or - indicate 
the DNA strand on which the region is encoded). Grey bars in the 
background indicate the two-fold cut-off level. 
(TIF) 

Figure S4 Pair-wise comparisons of expression in the IGEc/c 

area of mfsR mutants by microarray analysis, (A) P. putida UWG 1 
{ICEck-AmfiR, strain 4322) compared to P. putida UWGl (IGEc/c, 
strain 2737), (B) P. putida UWGl {IGEdc-AmfsR-A'marR, strain 
3453) versus P. putida UWGl (IGEc/c, strain 2737). (G) P. putida 
UWGl {ICEclc-mfsR::Km^, strain 2961) versus P. putida UWGl 
(ICEc/c, strain 2737). Panels indicate comparisons of exponentially 
growing or stationary phase ceUs, with hybridization signals on the 
plus- (open symbols) or minus-strand (closed symbols) of IGEc/c. 
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Dots indicate the ^log-fold change of hybridization signal per 
microarray probe in the comparison, plotted at their distance 
along the ICEcfc sequence pC-axis; in kb). A scheme of ICEcfc is 
redrawn at the bottom of each section, with regions of interest as 
grey boxes (+ or - indicate the DNA strand on which the region is 
encoded). Grey bars in the background indicate the two-fold cut- 
off level. 
(TIF) 

Figure S5 Observed and expected population growth in cultures 
of P. putida with wild-type ICEcfc and the AmfsR mutant. (A) 
Measured turbidities in triplicate batch cultures growing on 3CBA 
as sole carbon and energy source. Note the increased lag time of the 
P. putida AmfsR mutant. (B) Modeled population growth of all cells 
(brown) and tc cells (red) in two subsequent batch cultures growing 
on 3CBA; the second being inoculated with 1 : 1000 volume from the 
first culture in stationary phase. Scc-narios show predicted behavior 
for wild-type ICEcfe (with a probability P,c of 0.025 of tc cell 
appearance in stationary phase) and for AmfsR (with Pf^ = 0.5). Note 
how a tc population appears in stationary phase, which is 
transferred to a new culture, but rapidly dies as a result of activation 
of the par-shi system [30] . Note further how exponential growth rates 
remain the same for both wild-type ICEt/c and the AmfsR mutant. 
Modeling based on Monod kinetics, using the following parameters: 
total C = 4.25 mg C/ml, dry weight of a single cell = 2' 10''^ g, 
yield = 0.3 g/g. Kg = 0.02 mg/mL, generation time of non-tc 
cells = 1 h, generation time of tc-cells = 2 h, death rate of tc 
cells = 0.4 h-'. 
(TIF) 

Figure S6 Population sizes of P. putida UWC 1 carrying wild-type 
(^Til , 2738) or mutant ICEc/c, measured as colony forming units 
(CFU) on 5 mM 3-chlorobenzoate agar plates, per ml of 
resuspended culture spotted and incubated for 48 h on 0.5 mM 
3-chlorobenzoate containing agar medium. (A) P. putida UWCl 
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